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Summary. Modeling of radio and X-ray observations of supernovae interacting 
with their circumstellar media are discussed, with special application to SN 1993J 
and SN 2002ap. We emphasize the importance of including all relevant physical 
mechanisms, especially for the modeling of the radio light curves. The different 
conclusions for the absorption mechanism (free-free or synchrotron self- absorption), 
as well as departures from an p oc r~ 2 CSM, as inferred by some authors, are 
discussed in detail. We conclude that the evidence for a variation in the mass loss 
rate with time is very weak. The results regarding the efficiencies of magnetic field 
generation and relativistic particle acceleration are summarized. 



1 Introduction 

The interaction of supernovae (SNe) with their circumstellar medium (CSM) 
offer important clues to both the nature of the SN progenitors, the hydrody- 
namics of the explosion, the environment of the SN, and the physics of high 
velocity shock waves. The now convincing connection of Type Ic SNe and 
GRB's has also made the study of the SN environment especially interest- 
ing. Basically, the standard picture of the SN interaction with the surround- 
ings is only a non- relativistic version of the standard afterglow scenario (e.g., 
Chevalier, this Proceedings). We will discuss a few issues related to the SN 
interaction. For a more detailed exposition see the recent review in jS] . 



2 The standard model 



Supernova progenitors come in basically two flavors; extended red supergiants, 
or compact, hot stripped stellar cores. The former are thought to be related 
to Type IIP SNe, while the latter are most likely related to Type lb and Ic 
SNe. Type IIL, Iln and lib SNe probably represent a decreasingly massive 
hydrogen envelope. It is also tempting to identify this as a sequence of in- 
creasing ZAMS mass, with the Type IIP representing the most common low 
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mass progenitors, and the Type Ic's originating from the most massive Wolf- 
Rayet stars. A complication is the fact that a large fraction of stars are in 
close binary systems. Mass transfer between the companions can in this case 
lead to strong mass loss even for stars of comparatively low mass. 

Massive stars have in general strong stellar winds. On the main sequence 
the blue supergiants have fast winds with a wind velocity of 500—3000 km s _1 . 
As the star evolves to the red supergiant stage the escape velocity decreases 
dramatically and consequently the wind velocity to 10 — 30 km s _1 . In most 
of this stage the mass loss rate is 10 -6 — 10~ 5 M Q yr _1 . There are, however, 
indications, both from the observations of SNe (see below), and from stellar 
evolution calculations 14 that the star in some cases can undergo a stage 
with a superwind, similar to what occurs in AGB's, with mass loss rates of 
~ 10~ 4 M Q yr^ 1 . The duration of this stage must be very short, ~ 10 4 yrs. 
Finally, if mass loss is important enough, the star may evolve to the Wolf- 
Rayet (WR) stage, with mass loss rate ~ 10~ 5 M Q yr _1 and a wind velocity 
of 1000 — 5000 km s _1 , depending on the evolutionary stage of the WR star 
|3()| . For single stars this occurs only for stars more massive than ~ 22 — 40 
Mq for solar metallicity, depending on rotation [2Z|- In a binary system this 
may, however, occur at considerably lower mass. Because the CS density is 
p cs = M /(47r-u w r 2 ), where M is the mass loss rate, u w the wind velocity and 
r the distance from the star, the CS density into which the SN explodes can 
differ by several orders of magnitude, depending on the evolutionary stage in 
which it explodes. 

The collision of the supernova ejecta with the surrounding gas generates 
a strong shock wave, which expands with a velocity 20 — 30% larger than the 
maximum velocity of the ejecta. The temperature behind this forward shock 
is ~ 10 9 K. The pressure behind the shock will send a reverse shock back into 
the ejecta. Because of the higher density, the velocity of the reverse shock will 
only be V cs /(n — 2) ^ 500 — 1000 km s~ x , depending on the density gradient of 
the ejecta, p C j cx r~ n [5]. For polytropic envelopes the outer parts of the ejecta 
have n ~ 10 |24j . Departures from such a structure may, however, lead to both 
steeper and shallower gradients. For typical ejecta gradients the temperature 
will only be 10 7 — 10 8 K behind this shock. Therefore, and because of the 
high ejecta density, cooling will in general be important for the reverse shock 
|3 E3j • This will create a thin shell of cool gas between the reverse and CS 
shocks, often referred to as the cool, dense shell. The forward shock will in 
contrast usually be adiabatic, unless the CS density is very high. 

During the last two decades radio, optical and X-ray observations have 
greatly added to our understanding of the SN environment (see the papers by 
Immler, Filippenko and Van Dyk). Of these, the radio provides the cleanest 
signature of CS interaction. Because of this, and also because the interpre- 
tation of these observations have generated some confusion, we will discuss 
these observations in some detail, as well as the complications going into the 
analysis. 
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2.1 Radio emission and absorption 

The radio emission arises as a result of relativistic electrons accelerated in the 
shock region and emitting synchrotron radiation in the magnetic field. The 
exact mechanism of the acceleration is not well known, although a first or- 
der Fermi mechanism would be natural across the shock. Particle acceleration 
may, however, also occur behind the shock, close to the contact discontinu- 
ity, separating the shocked ejecta and shocked CSM. The generation of the 
magnetic field is not much better understood, although there has been recent 
progress in this area. In particular, numerical simulations [2H] based on 
ideas of Medvedev & Loeb .26 for GRB's, have shown that the Weibel in- 
stability, which is a particular version of the two-stream instability, can give 
rise to a strong magnetic field at the shock. Although the simulations have 
mainly been done in the context of relativistic shocks, the same mechanism 
should work also for non-relativistic shocks. An important issue for this type of 
simulations is to determine the efficiency of conversion of the thermal energy 
of the shock into magnetic field and relativistic electron energy. The latter 
two are usually characterized by e c = it c /utherm and ee = UB/uthcrm, where 
uthorm = Q/&pV 2 . Most likely, these parameters depend on both the shock 
velocity, the CSM density and chemical composition. 

The radio emission may be affected by free-free absorption from the 
surrounding CSM and synchrotron-self absorption (SSA) by the same rel- 
ativistic electrons that emit the radiation. Assuming the relativistic elec- 
trons to be injected behind the forward shock with an energy distribution 
given by dn c ("yj/dy oc 7 _p for 7 m ; n < 7 < 7max, the synchrotron spec- 
trum is given by F v — nR 2 S u [l — exp(— t„)] where S v oc v 5 / 2 /B 1 '' 2 and 
t„ oc B 3 / 2+a N c v~ a ~ 5 / 2 . Here a = (p — l)/2 and N e is the column density of 
relativistic electrons. In the case of negligible cooling N c oc Rn e oc Ru c . At 
low frequencies the optically thick spectrum is given by 

F v oc R 2 u b ' 2 B- l /\ (1) 

and at high frequencies 

F v oc R 2 N c B 1+a v- a oc R 3 n c B 1+a iy- a . (2) 

A fit to the spectrum, covering the peak frequency (y = ^ pC ak), therefore 
allows a determination of both the magnetic field and the density of relativistic 
electrons, if the radius of the emitting region is known. As we discuss below, 
this can usually be determined from optical line widths or in rare cases from 
VLBI. 

From these expressions we can derive an estimate of the brightness tem- 
perature at the peak of the spectrum (i.e., where r ~ 1) given by 

T b « 8.2 x xl0 10 .g(p) 2 / 17 (£ c / eB ) 2/17 (^/10 29 erg Hz- 1 ) 1 ^, (3) 
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independent of the value of v pe &k- The value of g(p) varies slowly with the 
parameters specifying the energy distribution of the relativistic electrons, for 
example, g(2) = l/hi(7max/7min)- The corresponding Lorentz factor is 

7 P cak « 1.0 x 10 2 (e o /e B ) 2 / 17 (F„/W 29 erg Hz" 1 ) 1 / 17 , (4) 

showing that 7 pe ak is fairly insensitive to various parameters, and is expected 
to vary by less than an order of magnitude. 

The second possibility, free-free absorption, is decoupled from the emission 
region, and only depends on the properties of the CSM. For an p oc r~ 2 CSM 
density, the free-free optical depth is 

r ff oc v~ 2 (M/u w ) 2 T- 3 / 2 R- 3 . (5) 

Here, R a = V s t. To estimate the relative importance of these mechanisms we 
determine the frequency of optical depth unity. For free-free absorption, we 
obtain 

W oc (M/u^T^V^h- 3 / 2 . (6) 

With 

B 2 /8tt oc e B pV 2 oc e B (M / Uvj )V 2 / R 2 oc e B (M/u w )/t 2 (7) 
and N e oc e e R s B 2 / e B we obtain for the SSA frequency 

^ssa oc (e B ec) 1/3 (A// Ww ) 2 / 3 F s 1 / 3 t- 1 . (8) 

The ratio of these is 

WssA oc (e B e e )- 1 /3(M/ Uw )V3y-li/3 T -3/4 t -i/2_ (9) 

Therefore, the relative importance of the two absorption mechanisms de- 
pends strongly on the wind velocity of the progenitor, the efficiencies of pro- 
ducing relativistic electrons and the strength of the magnetic field, the CSM 
temperature and especially the shock velocity. The latter point has been em- 
phasized by Chevalier [§], who find that most Type Ic SNe fall into the SSA 
category, because of their high expansion velocities and high wind velocities. 

For a consistent modeling it is crucial to include all relevant energy loss 
mechanisms for the relativistic electrons. These include the usual synchrotron 
losses, Compton losses on the photospheric and/or the synchrotron photons, 
and Coulomb losses. While the former two mainly act on the high energy elec- 
trons, and steepen the spectrum, the latter is most important at low energies, 
and lead to a flattening of the spectrum. Even with a constant value for p, this 
will lead to a curved optically thin spectrum. In practice, the kinetic equation 
for the electron distribution should be solved at each energy. More details on 
these processes can be found in 
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When the cooling time starts to become comparable to the dynamical 
time (i.e., f), radiation losses will affect the emitted synchrotron spectrum. 
For synchrotron cooling 

tsynch/t = 0.47(t/10days)( 7 e B )" 1 (M/10- 5 M Q jO^K/IO kms^ 1 ). 

(10) 

Since 7 > f0 2 , synchrotron cooling is expected to be important for super- 
novae with a red supergiant progenitor, unless £b *C 1. However, cooling in 
supernovae with a Wolf-Rayet star progenitor cannot be excluded, since the 
cooling time may be shortened due to inverse Compton losses. This can occur 
in situations where e c 3> £b (synchrotron self-Compton) or for very luminous 
supernovae for which scattering of photospheric photons may become impor- 
tant. As mentioned above, knowing z/ pea k and R s makes it possible to deduce 
individual values for ee and e - The occurrence of cooling provides an extra 
constraint on the model; for example, the underlying assumption of a spheri- 
cally symmetric source geometry can be tested. The source properties can be 
further refined if the Compton scattered radiation is observed. Furthermore, 
if a low frequency flattening due to Coulomb losses is not properly account 
for, it is likely that model fitting will give a value of t'poak, and hence B, which 
is too large. 

If free-free absorption is important, the epoch of optical depth unity allows 
us to determine the important density parameter, M/u w . Because R s = V s t, 
one can determine the ratio M/u w from the time of optical depth unity at a 
given frequency, if the shock velocity and temperature of the CSM are known. 
The shock velocity can either be estimated from the maximum extent of the 
optical line widths, usually the Ha line, or in more rare cases, like for SN 
1993J, directly from the radius as inferred from VLBI observations. 

The temperature of the CSM is even more difficult to determine. At the 
time of shock break-out a burst of soft X-rays heats and ionizes the CSM. 
The spectrum and luminosity of this burst depends on the nature of the 
progenitor, and has recently been discussed by Matzner & McKee [21]. A 
compact progenitor has a very brief burst of hard radiation with a low total 
energy, while a red supergiant has a longer burst, but considerably softer and 
with a larger energy. The ring of SN 1987A offers a spectacular example of 
a circumstellar structure heated and ionized by the outburst of the SN. In 
this case the initial temperature of the ring was ~ 2 x 10 5 K, cooling on a 
time scale of years 122] ■ At the smaller distances relevant for the free-free 
absorption of the radio supernovae during the first months (10 15 — 10 16 cm), 
these effects will be even larger. 

Unfortunately, with the exception of SN 1987A, little direct information on 
the temperature of the CSM is available. Instead, one has to rely on theoretical 
calculations. Because of the low density the recombination and ionization time 
scales are long compared to the evolution of the ionizing spectrum and one 
has to calculate the evolution of the temperature and ionization of the CSM 
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from the out-break until the epoch of interest. This has only been done for 
a few cases, the Type IIL SN 1979C [20], SN 1987A [2T] |22] and SN 1993 J 
|12| . In all cases the CSM temperatures close to the shock are 10 5 — 10 6 K 
immediately after shock break-out. The temperature then decreases on a time 
scale of months, but can be > 2 x 10 4 K even a year after the explosion. This 
result is in sharp contrast to the simplified analyzes of the radio light curves 
which assume a constant temperature of in most cases ~ (1 — 3) x 10 4 K, 
and obviously leads to large errors in the estimated mass loss rates. Moreover, 
the decreasing temperature can mimic a change in the mass loss rate with 
time, as has been claimed for some of the radio supernovae. These effects are 
discussed in detail in |2"U] . 

3 SN 1993J 

3.1 Synchrotron self-absorption or free-free absorption? 

The radio observations of SN 1993J @] EHj are to date by far the best 
ones available. Equally important, the VLBI observations P EI] give direct 
information about the size of the radio emitting region, and therefore the shock 
velocity, as function of time. From this combination a detailed modeling of 
the spectra have been possible ^T] El CHI EE] • The conclusions with regard 
to the mechanism behind the absorption, as well as the structure of the CSM, 
however, differ substantially. 

The analysis in both J2] and [SSj were based on modeling the radio ob- 
servations by free-free absorption in a CSM. As was shown in ^T], this alone 
cannot reproduce the turn-off of the spectrum at low frequencies, which, es- 
pecially at late epochs, is considerably shallower than the exponential cutoff 
expected for free-free absorption. Instead, it agrees well with the standard 
F v oc v^l 2 spectrum expected from synchrotron self-absorption. Only at early 
epochs was free- free absorption found to be important, and then only for the 
longest wavelengths. 

The main reason why free-free is only marginally important is that the 
temperature of the CSM inferred from modeling of the heating and ionization 
by the radiation from the shock break-out is likely to have been very high, 
~ 10 5 — 10 6 K The parametric fits by Van Dyk et al., on the other 
hand do not determine the temperature either observationally or theoretically, 
and simply assume a constant temperature. Because rg oc T~ 3 / 2 , this partly 
explains the different importance of the free- free absorption in |SH] and |llj . 
Furthermore, Van Dyk et al. do not use the dynamical evolution from the 
VLBI observations. Finally, their fits do not include any cooling processes 
which affect the electron spectrum, and only assume a constant power law 
for this. Both Coulomb and synchrotron were shown to be important in JT] , 
which flatten the low energy spectrum and steepen it at high energies. The 
effect of this can be seen especially in the lack of agreement between the fits 
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Fig. 1. Magnetic field (left) and density of relativistic electrons (right) as a function 
of the shock radius for SN 1993J. The dashed lines show the expected evolution if 
the magnetic energy density and relativistic particle density scale with the thermal 
energy density, B 2 /8w oc pV 2 oc n rc i cx t~ 2 , while the dotted lines show the case 
when B oc r~ x and n re i oc r~ 2 |ll|. 



and the observations at the peak of the lowest frequency, caused by Coulomb 
cooling, which flattens the electron spectrum. 

The modeling in [2] is also based on SSA. They do, however, not include 
a self-consistent calculation of the electron spectrum, for example Coulomb 
cooling is not accounted for, and invoke an arbitrary cutoff of the electron 
spectrum at low energies. This explains the different quantitative conclusions 
from 1 1 1| : in particular, the neglect of Coulomb cooling in the model fitting 
causes the deduced value of the synchrotron self-absorption frequency to be 
too high. This, in turn, results in an overestimation of the strength of the 
magnetic field. 

A test of the SSA model was provided by the new low frequency observa- 
tions with VLA in 32 and the VLA and GMRT in 4 . The fluxes of both 
these sets of observations agreed well with the fluxes predicted in [TI]. In 
addition, the combined VLA and GMRT spectrum at 3200 days in 0] nicely 
showed the break in the spectrum caused by the synchrotron cooling to evolve 
as expected. 

One of the most interesting results of the modeling of SN f 993J was that 
not only could individual values for the energy densities in magnetic field 
and relativistic particles be derived but also their evolution with time could 
be determined (see Fig. 1). Both of these energy densities scaled with the 
thermal energy density behind the shock (i.e., ee and e c are constants) but 
the conditions are far from equipartition since cb ~ 0.14 and £ c ~5x 10~ 4 . 



3.2 X-rays 

X-ray emission from SN 1993 J were observed with OSSE on Compton/GRO 
[IE], ROSAT P3IH1, ASCA [SZj, Chandra [HH1 and XMM [Hj. During the first 
two months, the OSSE and ASCA observations showed a very hard spectrum 
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with kT ~ 100 keV. This agreed well with that expected from the forward, 
circumstellar shock - 1.1 x 10 2 (V;/10 4 km s" 1 ) 2 keV [P2] . 

Because of the proximity to the Sun, it could then not be observed until 
~ 200 days after explosion. The temperature was now only ~ 1 keV, and the 
column density had increased by a large factor [23135]. This transition was ex- 
plained in as a natural consequence of the radiative reverse shock and the 
presence of the cool, dense shell. As the SN expanded the column density of 
the CDS decreases cx t , and therefore gradually becomes optically thin. Be- 
cause the luminosity below ~ 10 keV is dominated by the reverse shock, there 
will be a hard to soft transition, once the CDS becomes transparent. From 
the X-ray flux the mass loss rate was estimated to be ~ 4 x 10~ 5 M Q yr -1 , 
in agreement with that estimated from the radio. 

3.3 Structure of the CSM 

There has been considerable confusion also with regard to the structure of the 
CSM of SN 1993J. In several papers [121 CHI EH1 EH| there have been claims 
of a CSM density varying as p cx r^ 1 - 5 ^ 1 - 7 . This has in turn been taken as 
evidence for a mass loss varying with time. Only in jll| was an p cx r -2 
found to reproduce the observations. It is therefore of interest to examine the 
arguments on which these quite different conclusions rest. 

As has already been discussed, the analysis in ^21 an d [HE] both neglected 
the effects of SSA, and are therefore not physically consistent. Their conclu- 
sions from the radio modeling with regard to the CSM should therefore be 
ignored. 

The X-ray analysis by Immler et al. j!5j assumes that the X-ray emission 
emerges from the forward, circumstellar shock. This is directly contradicted 
by the low temperature, ~ 1 keV, and high column density found from the 
ASCA observations later than 200 days [37] • As was discussed above, the 
X-ray flux and temperature at late phases are instead consistent with those 
expected from the reverse shock ^2 EH] ■ 

The analysis by Suzuki & Nomoto [3fi| is based on a hydrodynamical, 
consistent modeling of the interaction of the ejecta and CSM, and the results 
should therefore be taken seriously. The fact that they obtain a CSM density 
at small radii varying as p cx r~ 1,7 and a clumpy medium at large radii, 
however, depends on the ejecta structure they use. The specific model they 
use, 4H47, was designed to reproduce the early light curve, and had for this 
purpose to be mixed artificially. More detailed modeling ^B] also showed that 
these 1-D models are hydrodynamically unstable. A further problem with 
this model is that it does not reproduce the velocity evolution of the shock, 
as inferred by the VLBI observations. More detailed calculations of the X-ray 
evolution using an ejecta model which does reproduce the VLBI observations 
indeed show that a satisfactory reproduction of the X-ray observations can be 
obtained for an p cx r~ 2 CSM density. 
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In conclusion, the only self- consistent modeling of the CSM of SN 1993J 
are those in using radio observations and fAffi using X-rays. Because of 
the insensitivity of the radio observations to the ejecta structure, in contrast 
to the X-rays, we believe the former is the more reliable, and that an p oc r~ 2 
CSM density is the most likely. 

4 SN 2002ap 

Type Ic SNe are of special interest because of their relation to the GRBs. The 
recent SN 2002ap is together with SN 1998bw the best observed objects of 
this class, and we discuss some issues related to this. 

A special property of the Type Ic SNe is their high expansion velocities. 
Although this varies by a large factor from the relatively slow SN 19941 to 
SN 1998bw, they are in all cases much faster than the Type IPs and also SN 
1993J. A problem here is the fact that the expansion velocity of the radio 
emitting region is difficult to estimate directly from observations. The optical 
spectrum shows few clear line features at early time, implying that blending is 
very important. This only allows a lower limit to the velocity to be determined. 
For SN 1998bw this was - 60,000 km s _1 |2E|, and for SN 2002ap - 30,000 
kms" 1 !^. 

Because of the high expansion velocities SSA dominates the radio absorp- 
tion. The radio observations of SN 1998bw \T%\ have been discussed in |19j . 
The most important conclusions was the high expansion velocity, V B /c <; 0.9, 
and the large energy in relativistic ejecta. Unfortunately, the radio observa- 
tions of this and other Type Ic's did not allow a determination of the mass 
loss rate, because this depends on the unknown values of e and £b- While 
these are often assumed to be close to equipartition, there is little theoretical 
motivation for this. 

SN 2002ap is interesting in this respect because, in addition to good radio 
observations 0, it was also observed with XMM [23112]. The X-ray emission 
is most naturally explained as a result of inverse Compton scattering of the 
photosphcric radiation by the same relativistic electrons which are responsible 
for the radio emission [3j . This process also leads to cooling of the relativistic 
electrons, and thereby a steepening of the emitted spectrum, which agrees 
with the fact that the optically thin radio spectrum is as steep as F v oc i^~ - 9 . 

By modeling the radio light curves of the different frequencies one can 
determine the expansion velocity of the radio emitting region, as well as the 
relative values of e G and cb- We find that V s ~ 70,000 km s" 1 , and that e e 
and cb are roughly equal (i.e., the magnetic field and relativistic electrons are 
in approximate equipartition). The exact value of their ratio depends on the 
upper and lower cut-offs in the electron distribution, which are not directly 
observable. This is in contrast to SN 1993 J, where we found that ee 3> e c . It is 
also interesting to note that in both 1993J and 2002ap the injected distribution 
of relativistic electrons has p « 2, which is the theoretical value expected from 
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first order Fermi-acceleration at a non-relativistic shock. This is not directly 
apparent from the observations, since cooling alters the emitted spectrum so 
that the simple relation a = (p — l)/2 cannot be used to deduce the value 
of p. It could be that an electron distribution with p — 2 is valid for most 
supernovae and that the range of observed values for a is due to cooling. 

5 Conclusions 

The VLA radio observations of SN 1993J are unique in terms of both the tem- 
poral coverage and the quality of the observations. The fact that also VLBI 
observations exist adds to this characterization, and means that the size of the 
radio emitting plasma can be determined without further assumptions. This 
combination makes a detailed spectral analysis possible, and from this a de- 
termination of the magnetic field and relativistic electron density as functions 
of time. For this a consistent physical model is necessary, which includes both 
the radiative transfer and the effects of different energy loss processes on the 
electron spectrum. There is no need to invoke arbitrary parameterized models, 
which only hides the real physical parameters, and may result in misleading 
conclusions with regard to physical mechanisms as well as the structure of the 
CSM. 

Acknowledgement. We are grateful to Peter Lundqvist and Roger Chevalier for com- 
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